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1. Introduction. 

This paper is devoted to the Fokas-Lenells equation (FLE) [HE], [3], 

(1) 



u xy -\- u — 2iuvu x = 



v xy + v + 2iuvv x = 

which appears in physical applications in the form 

u xy + u =F 2i\u\ 2 u x = 0. (2) 

FLE is known (see pQ) to be an integrable equation describing the first negative 
flow of the integrable hierarchy associated with the derivative nonlinear Schrodinger 
(DNLS) equation. In papers [2j [3] the authors derived the FLE from the Maxwell 
equations for the propagation of femtosecond pulses in optical fibers, established its bi- 
Hamiltonian structure, elaborated the inverse scattering transform scheme for the case 
of zero boundary conditions and obtained the N-bright soliton solutions. The aim of 
this paper is two-fold. Firstly, I want to derive the N-dark soliton solutions. Secondly, 
I am going to establish the relationships between the FLE and other integrable models: 
the nonlinear Schrodinger (NLS) equation, the Merola-Ragnisco-Tu (MRT) equations 
[I] and the Ablowitz-Ladik (AL) model [5]. 

It this paper I present the commuting pair of the Backlund transformations (section 
[2]) which I use in section [3] to obtain the lattice representation of the FLE. After 
bilinearizing the lattice equations and establishing their relationship with the AL 
hierarchy I use the N-soliton solutions of the latter to obtain the N-soliton solutions 
for the FLE (section H]). Finally, in sectional discuss the links between the DNLS, 
NLS, MRT and AL equations. 

2. Backlund transformations. 

The method of this work is based on two transformations, T and 7", given by 

{v = —iv y + uv 2 
~ Vx V 

and 

{u = iu y + u 2 v 
X_ 1 (4) 
IU X U- 

These transformations possess two important properties: 

• Transformations T and T send solutions of the FLE to solutions of the FLE. 

• Transformations T and T are inverse when applied to solutions of the FLE. 

To prove the first statement consider the quantities A[u, v] and A[u, v] which are 
the left-hand sides of equations (CE|): 

A[u, v] = u xy + u- liuvux ^ 

A[u, v] = v xy + v + HuvVx 
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It can be shown by direct calculations that if pairs (u,v) and (u,v) are related by ([3]), 
then 

A[u,v] . d A [u, v] 

( d 2v\ 

A[u, v] = v 2 A[u, v] — I i— — I ) A[u,v}. (7) 



dy v J 

From these equations one can easily see that if u and v are solutions of the FLE, 
A[u, v] = A[u, v] = 0, then A[u, v] = A[u, v] = which means that u and v also satisfy 
the FLE. In a similar way one can show that A and A calculated for the T-transformed 
functions are given by 

A[u,v] = (i^- - — ) A[u,v] +u 2 A[u,v] (8) 
\dy u J 

Tr , A\u, v] . d A\u,v] 

A[ fi , fi ] = ^J + I --^i. (9) 

Thus the identities A[u, v] = A[u,v] = lead to A[u, v] = A[u, v] = 0. To summarize, 

f A\u,v] = A\u,v] = 

A[u,v}=Mu,v} = =► tL 1 n (10) 

I A[u, v\ — A[u,v\ = 

which proves the fact that transformations T and T send solutions of the FLE to 
solutions of the FLE. 
To demonstrate that 

fof =Tof= 1 mod (A, A) (11) 

consider the results of joint action of T- and T-transformations: 



u 



~ ~ / u \ \ v = —iv v + uv- 

*r TT U- \*-i-h (12) 



2 



and 

' u 




r , -rr[ ), { , t _ i i (is) 



IU X u 

where u, v, u and v are defined in ([3]), ((4]). Using definitions (Hj) one can obtain from 

u 1 

u — u = — A[u, v], v — v — — A[u, v] . (14) 

u x vv x u 2 x 

In a similar way definitions (Tl3|) and ([3]) lead to 

1 _ v 
u — u — — A[u, v], v — v = - - z Am, v\. (15) 

vi uv x u x 



These identities together with ( TTOl) prove ( TTTT) 

In what follows I will use the symbol T for both transformations, 

f = r, t = r- 1 (i6) 

bearing in mind that T-transformations are used only with solutions of the FLE. 
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3. Lattice representation of the FLE and bilinearization. 

Using the transformations discussed in the previous section I introduce (instead of single 
solution of the FLE, u and v) an infinite sequence of solutions u n and v n defined by 
iteration of T-transformations: 

^)=r-(:), „ = 0, ± 1,±2,... (17) 

In terms of u n and v n the definitions of transformations T and T -1 , equations (E]) and 
(@J, can be written as 

-id x v n - 



1 + UnVn-x 

id y v n = v n+ i- u n v 2 n 



and 

idrU 



U n +1 



u n+l v n {19) 

idyU n u n —i u n v n . 
It is easy to see that the above equations can be regrouped in two systems 



and 



idyV n = v n+ i - u n v z n 



id x ii n 
-id x v n 



U n +1 



1 + vlT Vn (21) 



1 + U n V n -i 

These equations are not new: they belong to the MRT hierarchy that has been 
introduced, in the framework of the reduction technique for Poisson-Nijenhuis structures, 
in [1] as an integrable lattice hierarchy, whose continuum limit is the AKNS hierarchy. 
In [1] the authors demonstrated that this hierarchy is endowed with a canonical Poisson 
structure and admits a vector generalisation. They solved the associated spectral 
problem and explicity contructed action-angle variables through the r-matrix approach. 
In this work I show that equations (1201) and (12TT) , as well as the whole MRT hierarchy are 
closely related to another discrete version of the AKNS equations, namely the famous 
AL hierarchy. To this end I bilinearize equations ( 1201) and ( 12T1) by introducing the 
r- functions p n , a n and r n related by 

T n ~ r n-lT n +l = PnO~n (22) 

as follows: 

u n = , v n = . (23) 

By straightforward algebra one ca obtain that 

id x u n - - Un+1 =^—W(+\-u n J n J (24) 

1 + U n+ iV n T n _iT n L J 
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id x v n 

1 + U n V n -i 

idyU n - M n _i + U 2 n V n 

-idyV n - v n+1 + u n v 2 n 



'n—X'n 
1 

1 

' n— 1 '■ n 



[2<-> - v n J n .i] 
r (+) _ r 

[K,^ - v n C n _i] 



where 



2<+) = iD x a n -r n - r n _icr„, + i 

X^ _) = iD x T n ■ p n - p n _iT n+ i 
Jn — iD x T n+ i ■ T n + p n <J n +l 

/C^ +) = iD y o n -r n - a n _ir„ +1 
KS n ^ = iDy r n • p n — r n _ip n+ i 
C n = iDy T n+ i ■ r n — o~ n p n+ i 

and D x and D y are the Hirota's bilinear operators defined as 
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(25) 
(26) 
(27) 

(28) 
(29) 
(30) 
(31) 
(32) 
(33) 



D\ D l y a{x,y)-b{x,y) 



Qk+l 



+ €,y + v)K x -£,y- v) 



d^ k drj l 

One can easily see that equations f )20|) and fl21|) can be satisfied by 
conditions 



.(34) 

t=n=0 

imposing the 



l n ±] = J n = = C n = 0. (35) 
Thus, instead of system ( 1201) and (12ip one can solve the bilinear system given by 



%D X 0~ n • T n 7~n—lO~n+l 


= 


(36) 


iD x r n ■ p n — p n ~lT n+ i 


= 


(37) 


iD x T n +i ■ T n + p n cr n +i 


= 


(38) 


zl^ cr n • r n — cr n -iT n+1 


= 


(39) 


iDy T n ■ p n — r ra „ip n+ i 


= 


(40) 




= 0. 


(41) 



Equations ( 136|) -( |4T|) belong to the AL hierarchy that has been introduced 
in terms of the functions 

_ &n _ Pn 

Qn j ' n 

equations (j3"6"|) - (j4"T|) together with f )22|) become 

id x q n = (1 - q n r n ) Qn+i 
-id x r n = (1 - q n r n ) r n _ x 



and 



id y q n = (1 - q n r n ) Qn-i 
-id y r n = (1 - g n r n )r n+1 



in [5]. Indeed, 
(42) 

(43) 

(44) 
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equations that describe the simplest flows of the AL hierarchy and that can be viewed 
as superposition of the famous discrete NLS and discrete modified KdV equations. 

That means that one can use the already known solutions of the AL hierarchy to 
construct the ones for the MRT equation and hence for the FLE. 



4. Dark soliton solutions for the FLE. 



In the appendix one can find dark-soliton solutions of the generalized AL equations [6], 
the system which is a generalization of (f36l)-fj4TT). Using these results solutions of the 
latter can be written as 

= T°, a n = T n \ Pn = (45) 

where 

C(x, y) = {is) m \ n2 /i m (x, y) v n (x, y) <(x, y) (46) 

and 

u™(x,y) = det\I + A(x,y)M m Af n \. (47) 
Here A(x, y) is a square matrix with the elements 

^*(*,w) = -r^ir, j,k = i,..,N (48) 

M -^ir^)- (49) 

L and R are diagonal matrices, L = diag (Lj) and R = diag (Rj), related by 

(L — I)(R — I) = s 2 I (50) 

(I is the A^ x A^-unit matrix) with parameters r and s being related by 

r 2 -s 2 = l. (51) 

The dependence of the solution on the variables x and y is given by 

a k (x,y) = 4° } exp {-ir (L k - R k )x - ir (L^ 1 - R k Y ) y) , (52) 

p(x,y) = fj,^ exp {— ir 2 (x + y)} 
v(x,y) = exp {— is 2 (x — y)} , 



(53) 



where ajj. , p^ and are arbitrary constants. 

Bearing in mind physical applications, one can impose on the matrices L and R 
some restrictions which lead to the involution 

v n = -rTn". (54) 

Note that this involution, which is 'natural' for the MRT equations, is different from 
a more usual one for the AL system (r n = — q^). However, the former is sufficient to 
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derive the 'physical' solutions for the FLE, v = —u, which is the main goal of the present 
paper. Thus, let us rewrite fH6|) imposing the restrictions 



r > 1 (55) 
(which implies reality of s), 

L k = r + sexp(i$ fc ), 

(56) 

R k = r + sexp(-^ fc ) , 

|z/(°)| = r \^°'\, and introducing the real matrix C, 

C = AM- 1 ' 2 . (57) 

This leads, after some trivial modifications of the above formulae, to the following result 
for the functions u n and v n given by fl23|) : 



u n {x, y) = is e-^V- 2 " f ' k \' y \ ——^ (58) 

det \5 jk + C jk (x } y)e< 2n+1 ^r k=1 



, s det \8 jk + C jk (x,y)e i{ - 2n ^- a ^\ N u . 

v n (x, y) = is e^r 2n l — — - ^ (59) 

det\8 jk + C jk (x,y)e^-^ k=1 

Here 8j k is the Kronecker's delta, 

<p(x,y) = x+ (2r 2 - l) y + p {0) , (60) 
cp^ is an arbitrary real constant, 



c 



(0) 



exp [2r (X k x - X k l y) sin f3 k ] 



C, t(x , y) = -* a (61) 



sm 

where c^. are arbitrary real constants, 

A fc = |r + se^ fe | (62) 

and 

a k = 2 arg (s + re l,?fe ) - 
= arg (r + se l,?fc ) . 

Comparing (I55|) and (1591 it is easy to see that they indeed satisfy (1511) . Finally, setting 
n = one comes to the dark-soliton solutions of the FLE: 



(63) 



. . . det \5,; k + C ik (x, y)e , , 
u(x, y) = is e-vW ,/ ^ - y =1 (64) 



det 

and 



\6 jk + C jk (x,y)e i f l '>\ j)k=1 



v(x,y) = -u(x,y). (65) 
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5. Intersections of hierarchies. 



To conclude this paper I would like to discuss the question of 'intersections' of different 
integrable hierarchies that one can find in the previous sections: the DNLS, AKNS, 
MRT and AL hierarchies. 

The FLE belongs, as has been shown in pQ, to the DNLS hierarchy, describing 
its first negative flow. Equations (1201) and (1211) belong to the MRT hierarchy, which, 
in its turn, can be viewed as a discretization of the AKNS hierarchy. Finally, the 
bilinearization of the FLE and MRT equations fl20l) and (|2ip has led us to the AL 
hierarchy. To summarize, in this paper we were dealing with four integrable models: 
DNLS, NLS, MRT and AL systems. Now I would like to present some explicit formulae 
expressing solutions of the first three equations/hierarchies in terms of solutions of the 
AL model, which is known to be rather 'universal' and to 'contain' inside many other 
integrable models (see lecture [7] or original papers [SI El [9]). To do this I need more 
equations from the AL hierarchy (its second positive flow) which can be presented in 
the bilinear form as 



D t On ■ T n 

A r n ■ p n 

Dt T~n+1 ' T n 



D„. <x, 



n+1 



D x T n+ i 

D x p 



ln-l 



Pn-1 



(66) 



a 



n+1- 



These equations are compatible with fl36|) - pT|) and fl22|) being the members of the same 
hierarchy. Thus the system consisting of fl36|) - pTT) . fl66l) and fl22|) has nontrivial solutions 
and one can prove by direct calculations the following results: 



For any n the functions 



Or. 



U=—, 



V 



Pn-1 



n 



constant 



7"n 1~n—l 

solve the DNLS equation 

iU t + U xx + 2iUV U x 
-iV t + V xx - 2iUV V x 
For any n the functions 

Q 



o 





R 



Pn-1 



n 



constant 



solve the NLS equation 

iQt + Q XX 
-iRt + R xx 
The functions 

0~ n -l 



2Q 2 R 
2QR 2 








u. 



Pn 

Tn-1 



(67) 



(68) 



(69) 



(70) 



(71) 



solve the MRT equations (j2D]) and Q2T}. 
These results can be extended to the level of hierarchies: 
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• For any n the functions U and V given by (loTl) where p n , a n and r n are solutions 
of the AL hierarchy solve all equations of the DNLS hierarchy. 

• For any n the functions Q and R given by ( 1691) where p n , a n and r n are solutions 
of the AL hierarchy solve all equations of the AKNS hierarchy. 

• The functions u n and v n given by (I7T1) where p n , <r n and r„ are solutions of the AL 
hierarchy solve all equations of the MRT hierarchy. 
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Appendix A. Dark solitons of the AL equations 



The dark solitons for the AL equations were obtained in [10] using the inverse scattering 
method. Here I give another, purely algebraic, derivation of these solutions which is 
based on the Fay-like identities for the determinants of some special matrices. 

To make the following formulae more readable and to expose the inner structure 
of equations (!36|) - (HT]) and (1221) I rewrite them in terms of the double infinite set of 
tau-functions, r™, given by 



' n ) 
.0 



and 



m+l 



-(m+l) 



T,, 



Pr, 



2 _ m m ] 



-(m-1) 



( T -n 



T n-l T n+l 



111 



1,2, 



(a.i; 



(A.2) 



In terms of the new tau-functions equations fl36|) - fj4T]) become 



iD x T™ +l 


'm 

■ T n = 


m m+l 
T n-l'n+l 


(A.3) 


iD x C +1 • 


T m = 

n 


m— 1 m+l 
n 'n+l 


(A.4) 


iDy 


■ T m = 

n 


m m+l 
'n+l'n-l 


(A.5) 


in r m ■ 


T m = 

n 


m— 1 m+l 
'n+l T n 


(A.6) 



subjected to the restriction (T22|) 



m— 1 m+l 



^m m 
'n— l'n+1- 



(A.7) 



The dark-soliton solutions of the AL hierarchy are made of the N x N matrices A 
that satisfy the equation 

LA — AR = \ £)(a\ (A.8) 

where L and R are constant diagonal matrices, | £ ) is constant A-column, | £ ) = 



(4, 



'AT J 



and (a | is A-row depending on the coordinates describing the AL flows, 
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(a(x,y)\ = (ai(x,y), ...,ajsf(x,y)). By straightforward algebra it can be shown that the 
determinants 

uj (A) — det |I + A\ (A.9) 
satisfy the Fay-like identity 

a (P - 7) + P(j — a) oj/3Uj ai + j(a - 0) uj^oj^ = (A. 10) 

where 

u = co (A) , w c = u (AfQ) , (AK^K V ) (AA1) 

with 



k c = (i-cl) (I- CRT 1 - 

In what follwos I use the doubly infinite family of matrices A™ defined by 
A^ +1 = A™M, A™ +1 = A™J\f 

where 

M = K r , U = K 00 
and relate the matrices L and R by 

(L - rl) (R - rl) = (r 2 - l) I 

which leads to 

KrKcn = Ki/ r . 

Now equation (1A.10I) with (a, 0, 7) = (r, 1/r, 00) can be rewritten as 

(. ,m\2 _ 2 m-l m+l 1 2 m m 

where 



and p is given by p = Vl — In a similar way ( 1A.10I) with (^,7) equal to (1/ 



(r, 1/r) and (r, 00) leads to the following identities: 

(ra - 1) < < +1 + < +1 £C = ra u%# 

(a - r) < +1 < + r(l - ra) < < +1 = ap 2 a;™" 1 

[a-r) u n+l u) n +r u) n u n+1 - a u n+1 u n . 

with 

u>Z = u (A™K a ) . 
Analysing (IA. 19[) ( p\.21[) with a = A(£) and a = where 



r m MO ~ 1] 



A.12) 
A.13) 
A. 14) 
A.15) 
A.16) 
A.17) 

A.18) 

r,oo), 

A.19) 
A.20) 
A.21) 

A.22) 
A.23) 
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and 

X(t]) — r 

v = M V ) [r\( v ) - 1] (A ' 24) 

(these equations play the role of the dispersion laws of the positive and negative 
AL subhierarchies under finite-density boundary conditions) one can conclude, after 
identifying the Miwa's shifts with the action of the K\^Ky X r and K^^K^ 1 , that it is 
possible to introduce the dependence of A™ on an infinite number of 'times' in such 
a way that become solutions of the functional equations closely related to the so- 
called functional representation of the AL hierarchy [5J. However, I do not discuss here 
these equations in details because my purpose is to solve only the simplest AL equations 
0A.3p -( fAT6]) which can be done as follows. 

Let the x-dependence of A™ being determined by 

id x A™ = A%X (A.25) 

where 



X = lim - 



= r(L-R). (A.26) 

Replacing in a ( lA.19j) - (lA.20j) with A(£) and taking the £ — > limit one comes to 

iD x • < = r 2 u™£ - uC (A.27) 

iD x < +1 • < = p 2 [< < +1 - ^ <+*] • (A.28) 
In a similar way, introducing the ^-dependence of A™ by 

id y AZ = AZY (A.29) 

with 

Y = lim - [K- m K; 1 -l]=r (L' 1 - R' 1 ) (A.30) 

?7->-o rj 



one can get, after taking the 77 — >■ limit in (1A.20j) -( lA.21j) with a = \{rf) 



iD y uC +1 • < = r 2 < +1 - < (A.31) 

uC + ! ■ < = P 2 K +1 "Ci 1 - < <fi] ■ (A-32) 
Now one can obtain solutions of the AL equations by simple transformation 

= p m2 r n2 ii m v n u™. (A.33) 
It is easy to show starting from (1A.17j) that p m2 r™ 2 factor ensures that tau-functions r, 



fit 

Tl 



satisfy (1A.7I) . Then, by simple calculations one can get that if p and v depend on x and 
V as 

p = pWexp{-ir 2 (y + x)} 

(Q) r • 2/ m I A. Ol J 

z/ = z/ w exp {— «p — x)\ 

then equations (TAlTri) -( lA~32l) imply that t™ solve AL equations (lA~3l) -( lAi)l) . 

Noting that determinants (1A.9|) are invariant under transformations A — > U^AU 
one can eliminate (without loss of generality the constants £j (£j — > 1) by redefining the 
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functions a k (x,y) (a k (x,y) — > a k (x,y)£k) coming to the expressions for the elements of 
the matrix A (jUJ). Also, considering solutions of the MRT equations and the FLE it is 
convenient to use instead of p another parameter, s, given by p = is (see section H]). 
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